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The generalized Fisher super-exchange antiferromagnetic model with uniaxial 
crystal-field anisotropy is exactly investigated using an extended mapping tech- 
nique. An exact relation between partition function of the studied system and 
that one of the standard zero-field spin-1/2 Ising model on the correspond- 
ing lattice is obtained applying the decoration-iteration transformation. Con- 
sequently, exact results for all physical quantities are derived for arbitrary spin 
values S of decorating atoms. Particular attention is paid to the investigation 
of the effect of crystal-field anisotropy and external longitudinal magnetic field 
on magnetic properties of the system under investigation. The most interest- 
ing numerical results for ground-state and finite-temperature phase diagrams, 
thermal dependences of the sublattice magnetization and other thermodynamic 
quantities are discussed. 
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1. Introduction 

The Ising model turned out to be very convenient model for the investigation 
of magnetic properties of strongly anisotropic materials. In particular, the two- 
dimensional spin-1/2 Ising model and other exactly solvable models play an im- 
portant role in statistical mechanics of cooperative phenomena. Although, in most 
of theoretical studies the exchange interaction is usually assumed to be of short- 
range, in many real magnetic insulators the exchange coupling between magnetic 
atoms is predominantly carried out by means of the indirect exchange through the 
intermediary non-magnetic atom (the so-called superexchange interaction) . In order 
to describe these materials theoretically. Fisher introduced a two-dimensional anti- 
ferromagnetic Ising model on the decorated square lattice, in which the interaction 
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between decorating magnetic atoms realized via intermediate non-magnetic atoms 
was considered [ d 121 ■ 

It is also well known that the Fisher superexchange antiferromagnet is one of 
the few exactly solvable models in the presence of an external magnetic field. In 
fact, until its establishing, the Ising model in a nonzero magnetic field was included 
in the class of unsolved problems in statistical mechanics. Consequently, the effect 
of external magnetic field on magnetic properties was investigated only applying 
various approximate methods, the solutions of which were in many cases in no good 
agreement with each other as well as with experiments. The typical example of the 
experimental as well as theoretical discussions is the behavior of initial susceptibility 
in antiferromagnetic solids. Some experiments and approximations predicted that 
susceptibility has a maximum at a critical point and falls away on either its site. 
According to others, the zero-field susceptibility passes through a maximum up to 
above the critical temperature. At the critical temperature it shows the infinite 
slope. 

Moreover, no approximate methods yield explicit expressions for the initial sus- 
ceptibility as a function of temperature, nor do they indicate how it behaves in a 
finite magnetic field. The only case in which the magnetic susceptibility as well as 
other thermodynamic quantities can be exactly calculated is Fisher two-dimensional 
model. However, this model is only one of a number of similar exactly solvable anti- 
ferromagnetic lattices and its simple generalization leads to a series of other models 
which could display a variety of interesting features. To the best of our knowledge, 
only two generalizations of Fisher model have been investigated so far, namely, the 
superexchange Ising model on the kagome lattice [E] and the decorated Ising model 
which is the combination of the Fisher model and the Shiozi-Miyazima model of 
dilute-Ising spin system [Hj. 

Considering this fact, the main purpose of the present work is to generalize the 
original superexchange model for arbitrary values of decorating atoms including also 
the effect of the uniaxial crystal-field anisotropy as well. 

The outline of this paper is as follows. In the next section the basic points of 
the exact solution are briefly presented. In Sec. 3 we discuss the most interesting re- 
sults for the ground-state phase diagram, critical temperature and thermodynamics 
properties of the system under investigation. Finally, some concluding remarks are 
given in Sec. 4. 

2. Model and method 

In this article we will consider a generalized Fisher superexchange antiferromag- 
net on the square lattice (see Fig.l), in which vertex sites are occupied by the fixed 
spin-1/2 atoms (open circles) and decorating sites by the spin-S* atoms {S > 1/2) 
(closed and obliquely lined circles). Hence, we assume that horizontal and verti- 
cal decorating spins interact with each other through the so-called superexchange 
leading to an antiferromagnetic long-range order. Similarly to Fisher [0112], we also 
suppose that all horizontal decorating spins are coupled with their nearest neighbors 
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through an antiferromagnetic interaction, whereas in the vertical direction, there is 
a ferromagnetic nearest-neighbor exchange interaction. 

Furthermore, including the effects of external longitudinal field and uniaxial 
crystal-field anisotropy on all decorating atoms, the total Hamiltonian of the system 
can be written in the form 

2N 2N 2N 2N 

n, = JY.stf^]-JT.sifil-DY.istf -HYst, (i) 

ig_B iS-B 

where and S^, respectively, de- 
note the standard spin operators of spin- 
1/2 and spin-S" Ising atoms, the first and 
second terms describe the exchange in- 
teraction between the nearest neighbors 
in the horizontal and vertical direction 
and the last two expressions describe the 
interaction of decorating atoms with the 
uniaxial crystal-field anisotropy D and 
the longitudinal magnetic field H, re- 
spectively. Finally, the parameter J > 
represents the exchange integral be- 
tween the nearest neighbors. 

For later convenience, it is useful to 
rewrite the total Hamiltonian in the 
form Hd = J2h=i'^h + E^i^t,, where 
the first (second) term represents the 
sum of bond Hamiltonians involving all 
interaction terms associated with the 
decorating atoms on the horizontal (vertical) bonds and that are defined as follows: 

Hk = aJSmi + i^l2)-SlH-D{Sl)' (2) 

with k = h and v for horizontal and vertical bonds, respectively. The parameter a 
specifies the type of exchange interaction, thus for k = h we put a = 1 (antiferro- 
magnetic exchange) and for k = v is a = —1 (ferromagnetic exchange). 

The most important point of our calculation is the evaluation of the partition 
function Zd of the investigated system. The validity of commutation relation for the 
bond Hamiltonians (i.e. [Hi,Hj] = 0, for i ^ j) enables one to rewrite Zd in the 
partially factorized form, namely 

N 

Zd = TYexp{-PHd) = Tr{^} J| ZhZ^. (3) 

{h,v) 

In the above, (3 = l/fcsT, (fc^ being Boltzmann constant and T the absolute tem- 
perature) and Tr{^} means a trace over the degrees of freedom of spin-1/2 Ising 
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Figure 1. Part of the decorated square 
lattice with vertex spin 1/2 (open cir- 
cles) and decorating spins 3/2 (closed 
and obliquely lined circles. J denotes the 
exchange interaction which is antiferro- 
magnetic (ferromagnetic) in the horizon- 
tal (vertical) direction 
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vertex spins. Finally, the bond partition functions and furnish traces over all 
remaining degrees of freedom and they are given by 

s 

Zfc = Tr^, exp(-/37i:fc) = exp{PDn^) cosh [pn[ J {fil, + fil^) - aH]) , (4) 

n=-S 

where Tr^^ denotes a trace over spin-S* decorating atom on the fcth horizontal {k = h) 
or vertical [k = v) bond, respectively. The latter relation implies the possibility of 
introducing the the decoration- iteration mapping transformation [ El El 

Tts, expi-pHk) = Aexp {pRfil^fil^ + (3HQk{fili + fil2)l^) k = h,v. (5) 

Considering all possible configurations of the spins fi^i and /i^2 previous equations 
one finds that Hq^ = —Hq^. Consequently, when equation (0) is substituted into 0, 
the magnetic contributions H^h and Hq^ belonging to the vertex spins of the lattice 
cancel out and the partition function of the system Zd reduces to the form 

Zd{(3,J,D,H) = A'^Zo{f3,R). (6) 

Here, Zq represents the partition function of the standard spin-1/2 Ising square lat- 
tice without external longitudinal magnetic field and the transformation parameters 
A and R are given by 

A = {W{J)W{-J)Wm'^' , ^^ = ^"{^^^W^}' 

where W{x) depends on the temperature, external magnetic field and spin S of 
decorating atoms and it is defined as 

s 

W{x) = exp{PDn^) cosh{pnx - pnH). (8) 

n=-S 

It is worth noticing that equation © relates the partition function of the studied 
model in the presence of an external magnetic field and that of the standard zero- 
field spin-1/2 Ising model on the square lattice. Since the exphcit expression for Zq 
is known [ 7J, we can straightforwardly calculate many relevant physical quantities 
based on the familiar thermodynamic relations. On the other hand, some important 
quantities, such as the staggered magnetization or correlation functions cannot be 
obtained within thermodynamic approach. Fortunately, this problem can be solved 
utilizing the following exact spin identities [IH] 

{f2{Sl, fill, fll^)) d 

where arbitrary function /i depends on vertex spin variables, the function /2 depends 
on the spin variables of the fcth bond only and the symbols (...)o and {--^d stand for 



(/i(A 



I 1 h^i 1 



' Trsj2(^fc, fill, Afc2)exp(-/??i;A 
, Tr5,exp(-/5?Tffc) 



(9) 

k = h,v. (10) 
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the standard canonical averages of the original and decorated lattice, respectively. 
For example, applying the identity (fTUj) . one simply attains the following results for 
the sublattice magnetization in the horizontal direction 

m'B,^ = {Sl)d = Aq + Aim\ + A2eA, (11) 

where coefficients Ai depend on the temperature, external magnetic field and crystal 
field. Moreover, we have introduced the sublattice magnetization m\ = {fil)d and 
correlation function ea = (AaAi^2)rf- Very similar equations but with different coeffi- 
cient can also be derived for the correlation function = {filS^d and quadrupolar 
momentum g| = {{S^Y)d with k = h,v. Although the derivation of final equations 
for the above-mentioned quantities is straightforward, the calculation by itself is 
lengthy and tedious. Therefore we do not present details here. 

Finally, let us briefiy comment on the choice of the order parameter in our 
model. In general, at low temperatures the antiferromagnets are distinguished by 
the long-range order with anti-parallel alignment of the nearest-neighboring mag- 
netic moments (or spins). This type of magnetic arrangement is usually described 
by the so-called staggered magnetization which is defined as a difference of reduced 
sublattice magnetization. It is also well-known that Ising-type antiferromagnets in 
higher dimensions undergo a second-order phase transition at some critical tempera- 
ture and the staggered magnetization is then nonzero below the critical temperature 
and vanishes above this temperature. 

However, as we have already mentioned at the beginning of this section, the 
long-range antiferromagnetic order between pairs of horizontal and vertical deco- 
rating atoms of the lattice is realized indirectly through intermediate non-magnetic 
vertex atoms. Moreover, the vertex atoms (that are not coupled to the external 
magnetic field) also exhibit a spontaneous magnetization below the transition tem- 
perature. For this purpose, it is convenient to define the order parameter of the 
model under investigation in the form of binary vector m = (|?72,^|, m\), where the 
first vector component represents the absolute value of staggered magnetization 
(m^ = [m^Bh ~ '^Bv)/'^) describing the long-range antiferromagnetic ordering of 
horizontal and vertical decorating atoms, while the second one describes the spon- 
taneous magnetization of the spin- 1/2 Ising sublattice. 

3. Numerical results and discussion 

Before discussing the most interesting numerical results, it is worth emphasizing 
that the superexchange model is generally slightly different from the standard anti- 
ferromagnetic Ising models on the square lattice usually treated by most authors. In 
particular, the standard antiferromagnetic model becomes ferromagnetic when the 
sign of exchange integral J is changed, whereas the superexchange model remains in- 
variant against the transformation J ^ — J. It is also worth noting that for 5=1/2 
and 5* = 1 our calculation recovers the results obtained by Fisher, and Mashiyama 
and Nara [01121 El- For this reason we concentrate here on the case 5 = 3/2 which 
has not been discussed yet. 
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Table 1. Values of sublattice magnetization, correlation functions and quadrupo- 
lar momentum for different ground-state phases. 
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Now, let us proceed with the discussion of the ground state. Possible phases 
at T = 0.0 can be identified throught the analysis of all sublattice magnetization, 
pair-correlation functions and quadrupolar momenta belonging to the decorating 
spins and those of original undecorated Ising lattice. The results of our investiga- 
tion are collected in table 1. As we can see, there exist six different phases, namely 
three superantiferromagnetic phases (SAPi, SAP2 and SAP), one field-induced fer- 
romagnetic phase (FIFP) and two paramagnetic phases (PPi and PP2)- For a better 
illustration, we have depicted in figure 2 the ground-state phase diagram in the D—H 
plane. 

Comparing the quantities evaluated 
above one easily observes that for suf- 
ficiently small magnetic fields, the dec- 
orating sublattice of the system al- 
ways exhibits an antiferromagnetic long- 
range order. This is characterized by 
the mutually opposite signs of sublat- 
tice magnetization belonging to them: 
m%j^ = —3/2 and m^^ = 3/2 in the 
SAPi, m%f^ = -1/2 and = 1/2 in 
the SAP2 and m^^ — —1/2 and m^^^ — 
3/2 in the SAP. Moreover, the sublattice 
magnetization m\ takes its saturated 
value in the whole region where the su- 
perantiferromagnetic phases are stable. 
Consequently, the spin- 1/2 Ising atoms 
localized in vertexes of the lattice exhibit here exclusively the perfect ferromagnetic 
arrangement. 

On the other hand, the magnetic field stronger than a certain critical value 
{He/ J > 1.0) overturns anti-parallel decorating spins to its direction and conse- 
quently destroys the long-range antiferromagnetic ordering in the system. Therefore, 
in the region of magnetic fields H > He, both sublattice magnetization m^^ and 
will be positive: n^%f^ — = 3/2 in the PPi, m^^ = — 1/2 in the PP2 
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Figure 2. Grond-state phase diagram 
in the D — H plane for the decorated 
superexchange antiferromagnetic model 
with S = 3/2. 
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Figure 3. Global phase diagram in the 
D — H — T space for the decorated 
superexchange antiferromagnetic model 
with S = 3/2. 



and m|j;, = 1/2, m|j^ = 3/2 in the FIFP. In fact, for H > the FIFP is the only 
stable ordered ground-state phase, which always exists in a relatively narrow region 
of the crystal fields. It is also interesting, that in spite of the external magnetic 

field being applied only on the decorat- 
ing atoms, its sufficiently strong value 
indirectly effects the alignment of the 
Ising spins in the vertexes of lattice. In 
fact, in both paramagnetic phases PPi 
and PP2 the vertex non-magnetic spins 
are strongly frustrated due to the com- 
petitive effect of the magnetic field, crys- 
tal field and antiferromagnetic exchange 
interaction. Indeed, the zero values of all 
relevant correlation functions (see table 
1) clearly confirm the presence of frus- 
tration. One also finds that at the exact 
phase boundaries, the relevant phases 
co-exist implying the possibility of the 
first-order phase transition at T = 0.0. 
Thus the ground-state behavior of our model is much more complex in comparison 
with the original Fisher model. 

Apart from the ground state properties, we have also calculated the critical tem- 
peratures at which the system undergoes the second-order phase transition. Our 
results for the global phase diagram in the D — H — T space are shown in fig- 
ure 3. As we can see, the critical temperature exhibits a very interesting behavior. 
Naturally, the ordered (disordered) phases are stable bellow (above) the critical tem- 
perature. Of course, in some regions the stability of disordered phases extends down 
to zero temperature in agreement with the ground-state phase diagram. One should 

also note that our model exhibits es- 
sentially the same critical behavior as 
Fisher and Mishiyama and Nara mod- 
els [ m El E! • This is obvious, since the 
decoration of the lattice by finite num- 
ber of atoms cannot lead to a change of 
the universality class of the model. As 
we have already mentioned above, one 
of the most interesting physical quanti- 
ties in antiferromagnetic models repre- 
sents the initial suceptibilty. In order to 
understand the behavior of this quan- 
tity, in figure 4 we have presented ther- 
mal variation of the initial susceptibility 
(i.e H = 0) for several typical values of 
the crystal-field anisotropy. Besides the 
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Figure 4. Temperature dependences of 
the reduced initial susceptibility for the 
decorated superexchange antiferromag- 
netic model with S = 3/2 and some typ- 
ical values of the crystal field. 
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Figure 5. Temperature dependences of the reduced components m%, m\ of the 

order-parameter and the reduced sublattice magnetization ^^/^^ the 
decorated supcrcxchange antiferromagnetic model with S = 3/2, D/J = —0.5 
and some typical values of the magnetic field. 



standard curves (see the curve labeled 
D/J — 0.5), we also obtained some new dependences with clear maxima in the 
low-temperature region (see curves labeled D/J = —0.6 and —0.38). Moreover, one 
finds that strong nagative value of the crystal field rapidly depresses the absolute 
value of the initial susceptibility and causes its very slow decreasing above Tc (see 
the case D/ J — —1.0). This kind of dependences cannot be observed in the original 
Fisher model, since they originate from the crystal-field effects. The curve obtained 
for D/J = —0.5 also represents a furter original result which exhibits a strong di- 
vergence as temperature approaches zero. This behavior indicates the appearence 
of the first-order phase transition, since in this case we are exactly located at the 
phase boundary of the ground state. 
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Figure 6. Temperature dependences of the reduced specific heat and susceptibility 
for the decorated superexchange antiferromagnetic model with S = 3/2, D/J = 
—0.5 and some typical values of the magnetic field. 



Finally, to illustrate the behavior of other physical quantities at finite tempera- 
tures, in figures 5 we have depicted thermal variations of the order parameter com- 
ponents m^, and magnetization m^^, m^^ of decorating spins. In both figures 
we have fixed D / J — —0.5 and have changed the values of magnetic field. Compar- 
ing the behavior of the relevant quantities in figure 5 one observes the characteristic 
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behavior of the order parameter components being nonzero below Tc and vanishing 
at all temperatures above Tc. On the other hand m^f^, exhibit nonzero values 
in the whole temperature region. Thus they cannot serve as good order parameters 
for the model. At T = 0.0 all quantities take the corresponding characteristic values 
depending on the strength of the magnetic field. Of particular interest is the curve 
for H/J = 1.0 which apparently indicates the co-existemce of the SAPi and FIFP 
at T = 0.0 in agreement with the ground-state phase diagram. The existence of 
different phases and competitive effect of the magnetic field, crystal-field anisotropy 
and temperature are also reflected in the thermal dependences of the specific heat 
and susceptibility that are shown in figure 6. Both quantities exhibit characteristic 
behavior with interesting low-temperature maxima and logarithmic divergence at 
the critical temperature. 

4. Conclusion 

In this work we have generalized the Fisher superexchange antiferromagnetic 

model for the case of arbitrary decorating spins. Excepting the external magnetic 
field wc have also investigated the effect of the uniaxial crystal-field anisotropy. Some 
preliminary results for the case of 5* = 3/2 decorating atoms have revealed very rich 
and interesting behavior. Further interesting results for the magnetization process, 
thermal and field dependences of the entropy, specific heat and susceptibility will 
be presented in the near future. Moreover, more realistic versions of this model 
including the transverse magnetic field or biaxial crystal-field anisotropy are also 
exactly solvable and may potentially enrich our knowledge about exactly solvable 
models in statistical mechanics. 
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